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ABSTRACT 

Type la supernovae (SNe la) play an important role in astrophysics, especially in the 
study of cosmic evolution. There are several progenitor models for SNe la proposed in 
the past years. In this paper, we have carried out a detailed study of the He star donor 
channel, in which a carbon-oxygen white dwarf (CO WD) accretes material from a He 
main sequence star or a He subgiant to increase its mass to the Chandrasekhar mass. 
Employing Eggleton's stellar evolution code with an optically thick wind assumption, 
and adopting the prescription of Kato & Hachisu (2004) for the mass accumulation 
efficiency of the He-shell flashes onto the WDs, we performed binary evolution calcu- 
lations for about 2600 close WD binary systems. According to these calculations, we 
mapped out the initial parameters for SNe la in the orbital period-secondary mass 
(logP 1 — M%) plane for various WD masses from this channel. The study shows that 
the He star donor channel is noteworthy for producing SNe la (i.e. ~ 1.2 x 10 _3 yr _1 
in the Galaxy), and that the progenitors from this channel may appear as supersoft 
X-ray sources. Importantly, this channel can explain SNe la with short delay times 
10 s yr), which is consistent with recent observational implications of young popu- 
lations of SN la progenitors. 
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1 INTRODUCTION 

Type la Supernovae (SNe la) are excellent cosmological dis- 
tance indicators due to their high luminosities and remark- 
able uniformity. They have been applied successfully in de- 
termining cosmological parameters (e.g. f2 and A; Riess et 
al. 1998; Perlmutter et al. 1999). It is widely believed that 
SNe la are thermonuclear explosions of carbon-oxygen white 
dwarfs (CO WDs) accreting matter from their companions 
(see a review of Nomoto et al. 1997). However, several key 
issues related to the nature of their progenitor systems and 
the physics of the explosion mechanisms are still not well 
understood (Hillebrandt & Niemeyer 2000; Ropke & Hille- 
brandt 2005; Wang et al. 2008; Podsiadlowski et al. 2008), 
which may raise doubts about the distance calibration being 
purely empirical and on the basis of the SN la sample of the 
low red-shift Universe. 

At present, two SN la explosion models are discussed 
frequently, that is, Chandrasekhar (Ch) mass model and 
sub- Chandrasekhar (sub-Ch) mass model. The synthetic 
spectra of the Ch mass model are in excellent agreement 
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with the early time spectra of most SNe la, while that of the 
sub-Ch mass model has difficulty in matching observations 
(Hoflich & Khokhlov 1996) . A CO WD can increase its mass 
to the Ch mass through a single degenerate (SD) scenario, 
where the CO WD accretes H/He-rich material from a non- 
degenerate companion star (Whelan & Iben 1973; Nomoto 
et al. 1984), or through a double degenerate (DD) scenario, 
where another CO WD mergers with it and the total mass 
of the two CO WDs is larger than the Ch mass limit (Iben & 
Tutukov 1984; Webbink 1984). Theoretically, it is suggested 
that the DD scenario likely leads to an accretion-induced 
collapse rather than a SN la (Nomoto & Iben 1985; Saio & 
Nomoto 1985; Timmes et al. 1994). 

For the SD Ch scenario, the companion is probably a 
MS star or a slightly evolved subgiant star (WD + MS chan- 
nel), or a redgiant star (WD + RG channel) (Hachisu et al. 
1996, 1999a, 1999b; Li & van den Heuvel 1997; Langer et 
al. 2000; Han & Podsiadlowski 2004, 2006; Chen & Li 2007; 
Meng et al. 2008; Han 2008). Meanwhile, a CO WD may 
also accrete material from a He star companion to increase 
its mass to the Ch mass, which is known as the He star 
donor channel in this paper. The study of Iben & Tutukov 
(1994) showed that WD + He star systems can be formed 
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via binary evolutions. Yoon & Langer (2003) have carried 
out the evolution of a CO WD + He star system with a 
1.0 M Q CO WD and a 1.6 M© He star in a 0.124 day orbit. 
In this binary, the WD accretes He from the He star and 
grows in mass to the Ch mass. It is believed that WD + 
He star systems are generally originated from intermediate- 
mass binary systems. Thus, this channel may explain SNe 
la with short delay times (Mannucci et al. 2006; Aubourg et 
al. 2008). At present, however, there are no comprehensive 
investigations on this channel. 

The existence of WD + He star systems is supported 
by some observations. For example, KPD 1930+2752 is sug- 
gested to be such a system (Maxted et al. 2000; Geier et al. 
2007). Another example is from the discovery of a He nova, 
V455 Pup (Ashok & Banerjee 2003; Kato & Hachisu 2003). 
Recently, Kato et al. (2008) presented a free-free emission 
dominated light curve model of V445 Pup, based on an op- 
tically thick wind theory (Kato & Hachisu 1994; Hachisu et 
al. 1996). The light curve fitting in their study showed that 
the mass of the WD is lager than 1.35 M©, and half of the 
accreted matter remains on the WD, leading to the mass 
increase of the WD. Thus, V445 Pup is suggested to be a 
strong candidate of SN la progenitor (Kato et al. 2008). 

The purpose of this paper is to study the He star donor 
channel for SN la progenitors comprehensively and system- 
ically, and then to determine the valid parameter spaces for 
SNe la, which can be used in binary population synthesis 
(BPS) studies. In Section 2, we describe the numerical code 
for the binary evolution calculations and the grid of the bi- 
nary models. The binary evolutionary results is shown in 
Section 3. In Section 4, we estimate SN la birthrate of this 
channel. Finally, a discussion is given in Section 5, and a 
summary in Section 6. 



2 BINARY EVOLUTION CALCULATIONS 

In WD + He star systems, the He star fills its Roche lobe 
at He MS or He subgiant stage, and then the mass transfer 
begins. The He star transfers some of its material onto the 
surface of the WD, which increases the mass of the WD as a 
consequence. If the WD grows to 1.378 M©, we assume that 
the WD explodes as a SN la. Here we use Eggleton's stel- 
lar evolution code (Eggleton 1971, 1972, 1973) to calculate 
the binary evolutions of WD + He star systems. The code 
has been updated with the latest input physics over the last 
three decades (Han et al. 1994; Pols et al. 1995, 1998). Roche 
lobe overflow (RLOF) is treated within the code described 
by Han et al. (2000). We set the ratio of mixing length to 
local pressure scale height, a — l/H p , to be 2.0. The opacity 
tables are compiled by Chen & Tout (2007) from Iglesias & 
Rogers (1996) and Alexander & Ferguson (1994). In our cal- 
culations, the He star models are composed of He abundance 
Y — 0.98 and metallicity Z = 0.02, and all calculations for 
the He stars are carried out without enhanced mixing (i.e. 
the overshooting parameter, <5 ov , is taken to be zero; see 
Dewi et al. 2002). In addition, orbital angular momentum 
loss due to gravitational wave radiation (GWR) is included 
by adopting a standard formula presented by Landau & Lif- 
shitz (1971), 



where G, c, Mwd and M 2 are the gravitational constant, 
vacuum speed of light, the mass of the accreting WD and 
the mass of the companion He star, respectively. 

Instead of solving stellar structure equations of a WD, 
we use an optically thick wind model (Kato & Hachisu 1994; 
Hachisu et al. 1996) and adopt the prescription of Kato & 
Hachisu (2004) (KH04) for the mass accumulation efficiency 
of He-shell flashes onto the WD. If the mass transfer rate, 
\M 2 \, is above a critical rate, M cr , we assume that He burns 
steadily on the surface of the WD and that the He-rich mat- 
ter is converted into C and O at a rate M cr . The unprocessed 
matter is lost from the system, presumably in the form of an 
optically thick wind at a mass loss rate M w i n d = I-M2I — M cr . 
The critical mass transfer rate is 



M C1 



7.2 x 10~ 6 (Mwd /Mq - 0.6) M© yr _1 , 



(2) 



based on WD models computed with constant mass accre- 
tion rates (Nomoto 1982). 

Following assumptions are adopted when \M 2 | is smaller 
than M cr . (1) If IM2I is less than M cr but higher than 
the minimum accretion rate of stable He-shell burning, M st 
(KH04), it is assumed that the He-shell burning is stable and 
that there is no mass loss. (2) If \M 2 \ is less than M st but 
higher than the minimum accretion rate of weak He-shell 
flashes, M lo „ = 4.0 x 10~ 8 M©yr~ 1 (Woosley et al. 1986), 
He-shell flashes occur and a part of the envelope mass is as- 
sumed to be blown off from the surface of the WD. The mass 
growth rate of WDs in this case is linearly interpolated from 
a grid computed by KH04, where a wide range of WD mass 
and accretion rate were calculated in the He-shell flashes. 
(3) If I-M2I is lower than Mi ow , the He-shell flashes will be 
so strong that no mass can be accumulated onto the WD. 

We define the mass growth rate of the CO WD, Mco, 



M C = VHe\M 2 \, 



(3) 



where f)H C is the mass accumulation efficiency for He-shell 
burning. According to the assumptions above, the values of 
VHe are: 







M cr ^ |M 2 | Si M st , 
Mst > \M 2 \ > M low , 
|M 2 | < M low . 



(4) 



d In J G r _ 32G 3 M WD M 2 (M WD + M 2 ) 
dt ~~ 5c5 a 3 : 



(1) 



We incorporate the prescriptions above into Eggleton's 
stellar evolution code and follow the binary evolutions of 
WD + He star systems. The mass lost from these systems 
is assumed to take away specific orbital angular momentum 
of the accreting WD. We have calculated about 2600 WD 
+ He star systems, and obtained a large, dense model grid. 
The initial mass of the donor stars, M 2 , ranges from 0.85 M© 
to 3.1 M©; the initial mass of the CO WDs, M WD , is from 
0.865 M© to 1.20 Mq ; the initial orbital period of the binary 
systems, P', changes from the minimum value, at which a 
He zero-age MS (He ZAMS) star would fill its Roche lobe, 
to ~ 316 days, where the He star fills its Roche lobe at the 
end of the Hertzsprung gap. 
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Figure 1. A representative case of binary evolution calculations, in which the binary system is in the weak He-shell flashes phase at 
the moment of the SN explosion. In panel (a), the solid, dashed and dash-dotted curves show M2, A-fco an d A^WD varing with time, 
respectively. In panel (b), the evolutionary track of the donor star is shown as a solid curve and the evolution of orbital period is shown 
as a dash-dotted curve. Dotted vertical lines in both panels and asterisks in panel (b) indicate the position where the WD is expected 
to explode as a SN la. The initial binary parameters and the parameters at the moment of the SN la explosion are also given in these 
two panels. 
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Figure 2. Similar to Fig.l, but the binary system is in the stable Hc-shcll burning phase at the moment of the SN explosion. 
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Figure 3. Similar to Fig.l, but the binary system is in the optically thick wind phase at the moment of the SN explosion. 
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3 BINARY EVOLUTION RESULTS 

3.1 Typical binary evolution calculations 

In Figs 1-4, we present three representative cases in our bi- 
nary evolution calculations according to the condition of the 
binary system at the moment of SN explosion, and display 
one special case for producing SNe la. In panels (a) of these 
figures, we show the M 2 , Mrjo and Mwd varing with time, 
and the panels (b) are the evolutionary tracks of the He stars 
in the Hertzsprung-Russell diagrams, where the evolutions 
of the orbital periods are also shown. 

(i) Case 1 (see Fig. 1): The binary system is in the weak 
He-shell flash phase at the moment of the SN explosion. 
The binary shown in this case is (M 2 , M wd , l°g(P'/day)) 
= (1.35, 0.9, -1.20), where M%, M WD and P 1 are the initial 
mass of the He star and of the CO WD in solar masses, 
and the initial orbital period in days, respectively. In this 
case, the He star fills its Roche lobe after the exhaustion 
of central He (it now contains a CO core), and this results 
in Case BB mass transfer 0. The mass transfer rate |M 2 | 
exceeds M cr soon after the onset of RLOF, leading to a 
wind phase, where a part of the transferred mass is blown 
off in an optically thick wind, and the left is accumulated 
onto the WD. After about 5 x 10 4 yr, Af 2 | drops below M cr 
but still higher than M st . Thus, the optically thick wind 
stops and the He-shell burning is stable. With the continu- 
ous decreasing of \M%\, the system enters into a weak He- 
shell flash phase after about 5 x 10 4 yr. The WD always 
grows in mass until it explodes as a SN la in the weak He- 
shell flash phase. At the SN explosion moment, the mass of 
the donor star is Mf N = 0.8030 Mq and the orbital period 
log(P SN /day) = -1.1940. 

(ii) Case 2 (see Fig. 2): The binary system is in the 
stable He-shell burning phase at the moment of the SN ex- 
plosion. The binary in this case is (M 2 , M W di log(P'/day)) 
= (1.40, 1.00, -1.10). The binary evolves in a similar way 
to that of Case 1, but it is in stable He-shell burning stage 



1 We distinguish case BB (Roche lobe overflow occurs after He- 
core burning but before carbon ignition) from case BA (in which 
mass transfer is initiated during He-core burning) (see Dewi et 
al. 2002). 



when the CO WD reaches 1.378 Mq. The binary param- 
eters at Mwd = 1.378 M are M 2 SN = 1.0070 M and 
log(P SN /day) = -1.1296. 

(iii) Case 3 (see Fig. 3): The binary system is in the op- 
tically thick wind phase at the moment of the SN explosion. 
The binary in this case is (M 2 , M WD , log(Pyday)) = (1.80, 
1.10, -0.60). Since the initial mass of the WD is more mas- 
sive than that of the other two cases above, the WD may 
grow in mass to the Ch mass more easily, that is, it is still 
in the optically thick wind phase when Mwd = 1.378 Mq, 
at which M 2 SN = 1.2047 M® and log(P SN /day) = -0.6232. 

(iv) Case 4 (see Fig. 4): This case differs from the three 
above. Before the SN explosion, the binary system experi- 
ences two mass transfer phases. The binary in this case is 
(M|, M WD , log(P/day)) = (1.25, 1.10, -1.30). Due to the 
short initial orbital period of the system, angular momentum 
loss induced by GWR leads the He star to fill its Roche lobe 
before the exhaustion of central He, resulting in Case BA 
mass transfer. The mass transfer continues to proceed until 
the mass donor starts to shrink below its Roche lobe, termi- 
nating the mass transfer phase. After it has exhausted He in 
its core, the He star expands and fills its Roche lobe again 
when it evolves to the subgiant stage, leading to Case BB 
evolution. The following evolution of this system is similar to 
Case 1. When the accreting WD grows to the Ch mass, the 
binary is in the He-shell flash phase, and the binary param- 
eters are M 2 SN = 0.6983 M Q and log(P SN /day) = -1.1848. 

3.2 Initial parameters for SN la progenitors 

Figs 5-7 show the final outcomes of about 2600 binary evolu- 
tion calculations in the initial orbital period-secondary mass 
(logP 1 — A/ 2 ) plane, where the filled symbols are for those 
resulting in SN la explosions. The different cases described 
in Section 3.1 are plotted with different symbols, that is, 
the filled squares, circles and triangles denote that the WD 
explodes in the optically thick wind phase, in the stable He- 
shell burning phase and in the weak He-shell flash phase, 
respectively. Some WD + He star systems fail to produce 
SNe la because of He nova explosions (which prevents the 
WD growing in mass; the crosses in these figures) or dy- 
namically unstable mass transfer (resulting in a common 
envelope; open circles in these figures). Note that, the left 
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Figure 5. Final outcomes of the binary evolution calculations in the initial orbital period-secondary mass (logP 1 , Afi) plane of the CO 
WD + He star system for initial WD mass of 1.2 Mq, The filled symbols are for those resulting in SN la explosions, that is, the filled 
squares, circles and triangles denote that the WD explodes in the optically thick wind phase, in the stable He-shell burning phase and 
in the weak He-shell flash phase, respectively (see Cases 1-3 in Section 3.1). Crosses indicate systems that may experience He novae, 
preventing the WD from reaching 1.378 Mq, and open circles are those under dynamically unstable mass transfer. 




boundaries of the initial log P 1 — M\ plane in these figures 
are determined by the minimum value of log P 1 , at which 
the He ZAMS star would fill its Roche lobe. 

The contours of initial parameters for producing SNe 
la are also presented in these figures. The left boundaries 



of the contours in these figures (Figs 5-6 and panel c of 
Fig. 7) are set by the condition that RLOF starts when the 
secondary is on the He ZAMS 0, while systems beyond the 

2 Note that, the upper parts of the left boundaries are constrained 
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Figure 8. Regions in orbital period— secondary mass plane 
(logP r , Mj) for WD binaries (at the onset of the RLOF stage) 
that produce SNe la for initial WD masses of 0.865, 0.90, 1.0, 1.1 
and 1.2Mq. The region almost vanishes for M^-, =0.865 Mq. 



right boundary experience mass transfer at a very high rate 
due to the rapid expansion of the He stars in the subgiant 
stage and they lose too much mass via the optically thick 
wind, preventing the WDs increasing their masses to the 
Ch mass. The upper boundaries are also set mainly by a 
high mass transfer rate but owing to a large mass ratio. The 
lower boundaries are constrained by the facts that the mass 
transfer rate M2 should be high enough to ensure that the 
WD can grow in mass, and that the donor is sufficiently 
massive so that enough mass can be transferred onto the 
WD to reach the Ch mass. 

There is a time delay (up to about 10 6 yr) from the 
formation of most WD + He star systems to the onset of 
RLOFs, and this time delay varies with component masses 
and initial orbital periods. Fig. 8 shows the contours at the 
onset of RLOF for producing SNe la in the log P T — M\ plane 
for various WD masses (i.e. M^ D = 0.865, 0.90, 1.0, 1.1 and 



mainly by a high mass transfer rate because of orbit decay induced 
by GWR and large mass ratio, which leads to much of the masses 
lost from the systems in a way of the optically thick wind. 



1.2 Mq), where P r and M\ are the orbital period and the 
mass of the He companion star at the onset of RLOF, re- 
spectively. Note that, the enclosed region almost vanishes 
for M{v D = 0.865 Mq , which is then assumed to be the min- 
imum WD mass for producing SNe la from this channel. If 
the parameters of a CO WD + He star system at the on- 
set of RLOF are located in the contours, a SN la is then 
assumed to be produced. Thus, these contours can be ex- 
pediently used in BPS studies. One can send a request to 
wangbo@ynao.ac.cn for the data points of these contours 
and the interpolation FORTRAN code for these contours. 

Note, the contours for SNe la between Figs 5-7 and Fig. 
8 have some differences, that is, Figs 5-7 are for initial pa- 
rameters while Fig. 8 for those at the onset of RLOF. GWR 
may change the orbital period during this time, especially 
for binaries with short orbital periods (i.e. less than 1 day 
in this paper). For example, the left boundary of the con- 
tour for = 1-2 Mq in Fig. 8 has shorter periods than 
that in Figs 5-7, because angular momentum loss via GWR 
is more in short period binary systems than that in long 
period systems. 



4 BIRTHRATE OF SNE IA 

Adopting a prescription similar to that of Hachisu et al. 
(1999b) and based on Fig. 8 of this paper, we can roughly 
estimate Galactic SN la birthrate from the He star donor 
channel by using equation (1) of Iben & Tutukov (1984), 



v = 0.2 Ag 



dM 



AlogAyr \ 



(5) 



where Aq, A log A, Ma and Mb are the appropriate ranges 
of the initial mass ratio and the initial separation and the 
lower and upper limits of the primary mass for producing 
SNe la in solar masses, respectively. We give the details of 
the calculations in the following. 

To estimate the birthrate of SNe la, we divide the initial 



WD mass of Mw n into three intervals, 



0.9 - 1.0M, 



0i 

1. 0-1. 1M , and 1.1-1. 2Mq (see also Hachisu et al. 1999b). 
We ignore the range of M^,d = 0.865 — 0.9 Mq, since its 
birthrate is too small to contribute to the SN la birthrate 
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Table 1. The birthrates of SNe la for different WD mass intervals. 





A log A 


M A 


M B 


Aq 




(M e ) 




(Mq) 


(Mq) 




(lO-Syr" 1 ) 


0.9 - 1.0 


1.0 x 2/3 


5.60 


6.63 


0.02 


0.03 


1.0 - 1.1 


1.8 x 2/3 


6.63 


7.58 


0.18 


0.31 


1.1 - 1.2 


3.7 x 2/3 


7.58 


8.48 


0.34 


0.82 



as seen in Fig. 8. We use equation (11) of Yungelson et al. 
(1995) to estimate Ma and Mb, i.e. the final core mass (CO 
WD mass) versus the ZAMS mass relation, 

+<"»(*£)". w 

where Mwd and Mi are the mass of the final CO WD and 
of the ZAMS in solar masses, respectively. In addition, sim- 
ilar to Hachisu et al. (1999b), we also use an approximation 
relation to obtain A log A, i.e. 

A log .4 ^ |AlogP r , (7) 

where A log P r is taken from the SN la region in Fig. 8 
and the factor of 2/3 comes from the conversion between 
the period and the separation. Take WD mass interval of 
1.1 — 1.2 Mq as an example, we find that SN la explosions 
occur for the ranges of Mi, ; = 7.58-8.48 M , M Ho ,i = 0.95- 
3.0 M and A log A = 3.7 x 2/3, where Mi, ; is the initial 
mass of the primary at ZAMS, and Mne,i is the initial mass 
of the He star. Using Eggleton's stellar evolution code, we 
can obtain the initial ZAMS mass Mi = 5.60 Mq for the final 
He-core mass of 0.95 Mq (i.e. He star mass). Considering 
the initial mass of the secondary at ZAMS is not more than 
that of the primary at ZAMS, we constrain the range of M24 
from 5.60 Mq to Mi,;, where M24 is the initial mass of the 
secondary at ZAMS. Thus, we obtain iAvd, 1.1-1.2 ~ 0.82 x 
lO^yr -1 by substituting A log ,4 = 3.7 x 2/3, Aq = 1 - 
(5.60/8.48) = 0.34, M A = 7.58 and M B = 8.48 into equation 
(5). The SN la birthrates for different WD mass intervals 
are summarized in Table 1. Note, for different initial WD 
masses, the birthrate tendency of this channel is different 
from that of WD + MS channel, that is, massive WDs for 
producing SNe la have a high birthrate for the WD + He 
star channel. This is because that massive WDs easily have 
massive He stars as their companions. 

Finally, the sum of SN la birthrates for three intervals 
(0.9-1.0 M Q , 1.0-1.1 Mq, and 1.1—1.2 Mq) gives v ~ 1.2x 
10 — 3 yr _1 in the Galaxy, which is lower than that inferred 
observationally (i.e. 3 — 4 x 10~ 3 yr -1 ; van den Bergh & 
Tammann 1991; Cappellaro & Turatto 1997). 



5 DISCUSSION 

In our work, we have not considered the influence of rotation 
on the He accreting WDs. The calculations of Yoon et al. 
(2004) showed that, if rotation is taken into account, He 
burning is much less violent than that without rotating. This 
may significantly increase the accretion efficiency (?7He in this 
paper). Meanwhile, the maximum stable mass of a rotating 



WD may be above the standard Ch mass (i.e. the super- 
Ch mass model; Uenishi et al. 2003; Yoon & Langer 2005). 
However, in this paper we mainly focus on the Ch mass 
explosions of the accreting WDs. 

Yoon & Langer (2003) (YL03) only considered a case 
BB evolution of a CO WD + He star system which could 
produce a SN la. In this paper, we systematically studed 
the He star donor channel for the progenitors of SNe la (in- 
cluding case BA and case BB binary evolution) and showed 
the contours of SNe la for various initial WD masses. The 
binary studied in YL03 is located in the contours of this 
paper. However, there are some differences between our as- 
sumptions and theirs in mass accumulation efficiency. This 
is mainly because that the effect of wind mass loss in YL03's 
calculation is taken into account based an empirical formula 
of Wolf-Rayet star mass loss, while we used an optically thick 
wind model (Kato & Hachisu 1994; Hachisu et al. 1996) and 
adopted the prescription of KH04 for the mass accumulation 
efficiency of the He-shell flashes onto the WDs. 

The estimated birthrate of this channel is ~ 1.2 x 
10 -3 yr -1 , which is comparable to that of WD + MS chan- 
nel (this is considered as an important channel to SNe la; 
Han & Podsiadlowski 2004). Thus, the He star donor channel 
should not be ignored when we study the progenitors of SNe 
la, and this channel may increase the contribution of SD Ch 
scenario to SN la birthrate in theory. Note, the estimated 
birthrate may be higher than that in a real condition, since 
the long orbit period (i.e. > 1 day in Fig. 8) systems consid- 
ered to produce SNe la in equation (1) of Iben & Tutukov 
(1984) may not contribute to SNe la. Moreover, Umeda et 
al. (1999) concluded that the upper limit mass of CO cores 
born in binaries is about 1.07 Mq. If this value is adopted 
as the upper limit of the CO WD, the birthrate of SNe la 
from this channel will decrease to be ~ 0.2 x 10 _3 yr _1 in 
the Galaxy. 

Some precursory studies (Iben & Tutukov 1994; YL03) 
suggested that WD + He star systems may appear as su- 
persoft X-ray sources (SSSs) before SN la explosions. The 
initial orbital periods of WD + He star systems producing 
SNe la for this channel range from ~ 1 h to ~ 200 days, 
which may explain SSSs with various orbital periods, we 
can estimate the X-ray luminosity of a WD + He star sys- 
tem by Lx ~ £h c |M2|, where £h c = 6 x 10 17 ergg _1 is the 
energy generation rate owing to He burning. A WD + He 
star system has luminosity around 10 37 — 10 38 ergs _1 when 
He burning is stable on the surface of the WD, consistent 
with that of observed from SSSs (L x = 10 36 - 10 38 ergs -1 ; 
Kahabka & van den Heuvel 1997). Note, strong He II A4686 
lines are prominent in the luminous SSSs (see Kahabka & 
van den Heuvel 1997). Thus, we emphasize that SN la pro- 
genitors in the He star donor channel may appear as SSSs 
during the stable He-shell burning phase without winds 0. 
By multiplying the birthrate, ~ 1.2 x 10~ 3 yr -1 , with an 
average duration of the SSS phase, ~ several times 10 s yr, 
we estimate the current number of this type of progenitors 
as SSSs in the Galaxy to be from a few to several hundred. 
However, we do not observe such number of WD + He star 
systems in SSS phase. This may attribute to the Galactic in- 



3 Massive wind may absorb soft X-rays, which then cannot be 
detectable (Hachisu et al. 1999b; Kato & Hachisu 2003). 
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terstellar absorption of SSSs (Di Stefano & Rappaport 1994; 
Hachisu et al. 2008a). Moreover, Di Stefano & Rappaport 
(1994) also suggested that circumstellar matter may play 
a role in the obscuration of X-rays (see also Hachisu et al. 
2008a). 

Generally, a CO WD + He star system is produced from 
an intermediate-mass binary, which can be found in stellar 
populations with relatively recent star formation. The mini- 
mum He star mass for producing SNe la from this channel is 
0.95 Mq, with which we can roughly estimate the maximum 
delay time of this channel. The initial ZAMS mass Mi = 
5.6 Mq can produce He star mass of 0.95 M (see Sect. 4). 
The MS lifetime is about 81 Myr for a star of 5.6 M (Eggle- 
ton 2006). Moreover, we should consider the MS lifetime of 
He stars, since the He stars in long orbital period systems 
experience RLOF after the exhaustion of central He. The 
MS lifetime of 0.95 Mq He star is about 18 Myr. Thus, the 
maximum delay time of this channel is about 10 8 yr, i.e. the 
He star donor channel can explain SNe la with short delay 
times (< 10 8 yr; Mannucci et al. 2006; Aubourg et al. 2008). 
Note, Hachisu et al. (2008a) investigated new evolutionary 
models for SN la progenitors, introducing the mass-stripping 
effect on a MS or slightly evolved companion star by winds 
from a mass-accreting WD. The model can also explain the 
presence of very young (< 10 s yr) populations of SN la pro- 
genitors (see also Hachisu et al. 2008b). 

When WD + He star systems explode as SNe la, the SN 
la spectra may more or less show sign of He. Unfortunately, 
He associated with SNe la fail to be detected in the past 
years (Marion et al. 2003; Mattila et al. 2005). Therefore, 
Kato et al. (2008) suggested that SNe la from this chan- 
nel may be rare. However, further study on this channel is 
necessary, since this channel can explain some SNe la with 
short delay times. In addition, GWR is strong in WD + 
He star systems with short orbital periods. Thus, the WD 
+ He star systems with short orbital periods are possibly 
detectable sources of GWR. 



6 SUMMARY 

Using an optically thick wind model (Kato & Hachisu 1994; 
Hachisu et al. 1996) and adopting the prescription of KH04 
for the mass accumulation efficiency of the He-shell flashes 
onto the WDs, we systematically studied the He star donor 
channel for producing SN la progenitors and showed the 
parameter spaces for SNe la from this channel. The study 
shows that this channel is noteworthy for producing SNe 
la, that is, the birthrate of SNe la is ~ 1.2 x 10~ 3 yr -1 
in the Galaxy. In addition, the minimum mass of CO WD 
for producing SNe la from this channel may be as low as 
0.865 Mq. According to the orbital period range and X-ray 
luminosities of WD + He star systems, we considered that 
WD + He star systems may appear as SSSs before SN la 
explosions, which is consistent with YL03. Importantly, the 
He star donor channel can explain SNe la with short delay 
times (< 10 8 yr). Meanwhile, WD + He star systems with 
short orbital periods may be detectable sources of GWR, 
and the contours of log P T — M| plane at the beginning of 
RLOF can be used in BPS studies. In future investigations, 
we will study the properties of companions at the moment 



of SN explosion, which might be verified by future observa- 
tions. 
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